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ABSTRACT 

We calculate the one-loop QCD corrections to t — > tiXj using dimensional reduc- 
tion scheme, including QCD and supersymmetric QCD corrections. The analytic 
expressions for the corrections to the decay width are given, which can easily be 
extended to t — > Xjbi- The numerical results show that the correction amounts to 
more than a 10% reduction in the partial width relative to the tree level result. We 
also compare the corrections in the no-mixing stop case with those in the mixing 
stop case. 
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1. Introduction 



The top quark has been dicovered by the CDF and DO Collaborations at the Fermilab 
Tevatron[l]. In the Standard Model t — > W + + b is the dominante decay mode. Beyond 
the SM, in addition to the top decay into possible charged Higgs bosons plus bottom, a 
potentially important decay channel of the top quark is the supersymmetric decay into a 
lighter stop plus a neutralino, which has been extensively discussed at tree level [2]. It is 
generally expected that the lighter of the two stops is significantly lighter than the other 
squarks because the large top quark Yukawa coupling drives the diagonal stop masses to 
small values and enhances the off-diagonal mixing of left-handed and right-handed stops, so 
the present squark mass limits would not apply to the lighter stop. The best current lower 
bound on the stop mass is 55GeV and comes from LEP, operating at yt~s = 130 — 140 GeV[3]. 
The DO experiments at the Tevatron have excluded the existence of a stop lighter than 
100 GeV, albeit under certain assumptions [4]. Since the lightest neutralino is the lightest 
supersymmetric particle the decay t — > t±Xi could occur in a reasonably large volume of 
the parameter space with a sizeable branching ratio [2]. The one- loop radiative corrections 
to both t — > W + b and t — > H + b have been calculated [5,6] but the radiative corrections 
to t — > tix'j and t — > Xjbi have so far not been calculated. In this paper we present the 
calculation of the one-loop 0(a s ) corrections to the top quark decay into the lightest stop 
plus a neutralino, including both QCD and supersymmetric QCD contributions. Our results 
can be generalized straightforwardly to the decay t — > Xjbi, where b\ is a light sbottom. 



2. Tree-level 

In order to make this paper self-contained we first present the relevant interaction La- 
grangians of the Minimal Supersymmetric Standard Model (MSSM) and the tree-level decay 
rates for t — > tiXj- The interactions of top and stop with neutralinos and gluinos are given 
by the Lagrangians[7] 

C m o = -y/it(LijP L + RijP^XjU + h.c, (1) 
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and 



where 



£tti9 = -9sT a t(ai - ba 5 )g a ti + /i.e., (2) 



1 1 

ai = -^(cos# — sin#), a 2 = Wcos# + sin#), (3) 

V 2 y2 

fex = — ^(cos# + sin#), 6 2 = — ^=(sin# — cos#), (4) 

V 2 V 2 

Lij = Aj cos — Cj sin 0, L 2 j = — Aj sin 9 — Cj cos 9, (5) 

= -Bj cos 9 — Dj sin 9, R 2 j = —Bj sin 9 — Dj cos (6) 



with 



A- - D* 9mtN ^ B -C* + 9 -^ (7) 



2 ATf * 2gS\ 



C 3 = ZeN^-litAN*, ( 8 ) 



and 



N' 3l = N 3l C w + N j2 S w , N' j2 = -N n S w + N, 2 C w , (9) 

Here S w = sin0 w ,C w = cos#w, Pl,r = |(1 =F 75), and iV^ are the elements of the 4x4 
matrix iV defined in Ref.[7], which can be calculated numerically. T a = \ a /2 are the Gell- 
Mann matrices and 9 is the mixing angle between left- and right-handed stops which are 
related to the mass eigenstates U in Eqs. (1) and (2) by 



/ ti \ _ ( cos 9 sin 9 \ ( t L \ 
\ t 2 J ~ \ - sin 9 cos 9 J \ t R J ' 

This rotation matrix, Eq. (10), diagonalizes the stop mass matrix[8] 



(10) 



M 2 = ( M I L + m 2 + 0.35cos(2/3)M 2 -m t (A t + //cot (3) \ 

1 \ -mt(At + n cot /3) Mf R +m 2 + 0A6cos{2(3)M 2 )' [ } 

where M~ l , M 2 r are the soft SUSY-breaking mass terms for left- and right-handed stops, /j, 
is the supersymmetric Higgs mass parameter in the superpotential, A t is the trilinear soft 
SUSY-breaking parameter, and tan/5 = v 2 /v 1 is the ratio of the vacuum expectation values 
of the two Higgs doublets. 
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The tree-level Feynman diagram for the decay t 
tree-level partial decay width is given by 

1 



tix°j is shown in Fig. 1(a), and the 



r n = 



,-A 1 / 2 ! 



2 2 2 ^ 



where X(x, y, z) 



+4Re(LljRij)m t m^q 
(x — y — z) 2 — Ayz. 



T I 2 



R ij\ 2 )i m t + m h -ml) 



12) 



3. Virtual corrections 

Since the conventional dimensional regularization violates supersymmetry, in our calcu- 
lation we will use dimensional reduction technique [9], which preserves supersymmetry, for 
regularization of the ultraviolet divergences in the virtual loop corrections, although there 
is only a small difference between the both schemes to first order in the QCD and weak 
couplings. In fact, in dimensional reduction scheme, at the one-loop level the e-scalars con- 
vert the dimensional regularization result to the result which would be obtained by simply 
performing the numerator algebra in four dimensions [9]. To regulate the infrared divergences 
associated with soft and collinear gluon emission we will give the gluon a small finite mass 
A which is legitimate for our purposes since the non-Abelian nature of QCD does not show 
up in this order. We will also adopt the on-shell renormalization scheme [10] in which the 
coupling constant and the physical masses are chosen to be the renormalized parameters. 
The finite parts of the counterterms are then fixed by the renormalization conditions that 
the quark and the squark propagators have poles at their physical masses. For the QCD 
and SUSY-QCD corrections to the decay t — > t±Xj, which we are considering, only the top 
quark mass and the stop mixing angle in the bare coupling need to be renormalized. By 
introducing appropriate counterterms the renormalized amplitude can be expressed as 

M ren = -iV2u(xj)(aP R + bP L )u(t) (13) 

with 

a = L^ + SL^ + Ll^Z^ + ^Z^ + Ll^ + A^ + A 5 ^-' 300 , (14) 

b = Rl j + SRl j + Rl j (^SZ t L + ^SZ 11 ) + R; j SZ u + A < i CD + A s L USY - QCD : (15) 
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where k Q L c R D and A S L U R S Y ~ QCD are the vertex corrections from the irreducible vertex diagrams, 
expressions for which will be given below. 5Ly and ■ are the shifts from the bare couplings 
to renormalized couplings and , as mentioned above, can be found by renormalizing the top 
quark mass and the stop mixing angle: 

*(m t ) $ m t 



8L\ 3 = L* 2j 56 + L ±j 



R,2jS9 + 



m t 



m t 



Y *(mt) 



-D*sm9. 



(16) 

(17) 
(18) 



A* cos9, i^ mt) 

The counterterms and the renormalization constants in Eqs.(14)-(17) are defined by 

m' = m t + 5m t , (19) 

e° = e + se, (20) 

t° = zl /2 t=(l + 5Z t L P L + 5Z t R P R ) 1 /H, (21) 

and 

t? = (1 + SZ u ) 1/2 ii + 5Z 12 t 2 . (22) 
Calculating the self-energy diagrams for the top quark in Figs. 1(b) and 2(a) we obtain 



m t 

SZ, L = 



u s Cf 
a s C F 

47T 



-2A + 4F { " 9) 



2A + 2Ff 3) +m 2 t {4G ( ; 



(ttg) m 9„, jpitgti) 



2F r> ~ -^auFZ 
m t 



i(tgti) 



(ttg) 



+2m 2 t a il G ( f gk) + 2m t m~ g a^l 



(tgti) 



(23) 



(24) 



and 



5Z« = 



47T 



2A + 2Fl ttg) + m 2 t (4G ( ; 



(ttg) 



(tgk 



+2m 2 t a u Gf 9k) + 2m t m- g a ll Gf h) 
where the sum over i(— 1, 2) is implied, and 





= didj + bibj, 




= didj-bibj, 


\j 


= a,jbj + djbi, 


i?(ijk) 
n 


= f dyy n log 
Jo 



™>iv(y - !) + -y) + m W 



(25) 

(26) 
(27) 
(28) 

(29) 



and 



Q{ijk) 



- I' dy- 

Jo ? 



y n+1 (i-y) 



(30) 



™iV(v - 1) + mj(l - y) + m\y' 
Here, A = - — 7 S + log47r with ^ E being the Euler constant and D = 4 — 2e is the space-time 
dimension. The color factor Cf = 4/3 for 577(3) and fi is the 't Hooft mass parameter in 
the dimensional regularization scheme. Similarly, from Fig. 1(c), 2(b) and 2(d), one finds, 
for the stop, 



5Z 



ii 



OL s Cl 
47T 



-F, 



(tihg) 



o 



-2F± 



^ - 2ml(G { ilil9) + Gf liig) ) 
m 2 t G { hm + m t m~ g sm(26)G { hm 



(31) 



and 



5Z 



12 



An 



cos(2^) 



4m t mg 
m~ — m~ 

(2 ti 



+ sin(2^) A 



2 2 

m- — m~ 



sin(2^)(A (m fi ) - A (m h )) + Am t m~ g F^ 



(tits) 



with 



A (m) 



(32) 



(33) 



We have fixed the wave function renormalization constants and the top quark mass coun- 
terterm by the on mass-shell renormalization scheme condition. The mixing angle countert- 
erm is fixed by the requirement that 58 exactly cancel the remainder of the sum of all the 
ultraviolet (UV) divergent terms in the square of the renormalized ampltude, insuring the 
UV finiteness of the physical observables. From this requirement we found that the mixing 
angle counterterm simply is the negative of the counterterm 5Zi 2 ; that is, 



56 = -5Z 



12- 



(34) 



This condition insures that all the ultraviolet divergences will cancel in the virtual corrections 
to the decay width, as will be seen below, and is in agreement with Ref . [1 1] . 

The calculations of the irreducible vertex corrections from Fig. 1(d) and 2(c) results in 



A QCD = A QCD Pl + a QCD Pr 

{(L* 1:j P R + Rl J P L )[A + 4C 24 



l L 1 L 

a s C F 



An 



+m 2 (2C + 2C n - C 12 + C 21 - C 23 ) + m? (2C + 2C n 

+C12 + C23) — m|o(2Co + 2Cn — C12 — C 22 + 6*23)] 

+2(L*.P L + ^^^^^(Cn - C 12 )} (-p, fcx, A,m t ,m fl ), (35) 



and 



jySUSY-QCD _ j^SUSY-QCD + ^SUSY-QCD 

= ^-^ { [{L 2j - L 2j cos 29 - L l3 sin 2fl)P L 

47T 

+(-P 2j - P 2j cos 26 - R Xj sin 2#)Pr] A 
+ [Sf [4C 24 + m 2 (C 21 - C 23 + C u - C 12 ) 
+m\(C 22 - C 23 ) + m|o(C 23 + C12)] 

+m t 2 Sj?(C 11 - Ci 2 + Co) + m t m i o[S , j? ) Ci 2 + ^(Q, + Cu)] 
+m^oSjJ ) (C + C 12 ) + m t m 5 [Sj? C + S$\C U - C 12 )]) P fi 
+ (4 2) [4C 24 + m 2 (C 21 - C 23 + C n - C 12 ) + m| (C 22 - C 23 ) 
+m|o(C 23 + C 12 )] + m 2 t Sf (C n - C 12 + C ) + m t m x o[S'jJ ) C 12 
+4 1} (C + C n ) + m^oS'jf (C + C 12 ) 

+m t mp[5 , j- ) C + SjJ^Cn - C12)]) Pl] (~P, fa, m iv mg, m t ), (36) 

respectively, where the sum over i(— 1,2) is implied. In Eqs.(35) and (36) 

SjJ — (oiii + (3ii)Rij, SjJ — (a>u — [3ii)Lij, 
Sj? = ( a u + Pu)Lij, Sj^ = (an — [3ii)Rij, 



Sjf — (fii — Aij)Ljj, sjf — (<7ij + \\i)Rij 
r( 7 ) — ^ . \ w c( 8 ) 



S'jj — (<7ij + Aij)Ljj, — (uij — \ii)Rij : (37) 

where = — frjOj, and Co, C^ are the three-point Feynman integrals given in the 
appendices of Ref. [12]. 

The virtual correction to the decay rate is then 

5T virt = \ 1/2 (m 2 t , m|o, m|)Pe {2(m 2 + m|o - m?) [(L y L^ + RijR* 2j )(S9 + 5Z 12 ) 



+^(\L lj \ 2 5Z t R + \Ru\ 2 5Z t L ) + {L l3 Sf ] + i^-Sf Vo^ Y ~ QCD + ^iA + R l3 8 2 



J 

*(m t ) , p T *(mt)\Srrit 



(L l3 R* 2j + R^L^Se + 5Z 12 ) + (Lyi^ + R l3 L^ 



+ (Li 3 Rl, + R^qSZn + 8% CD ) + -(L^SZf + RuLtfZ?) 



Uz ll + 5^ D ) + \ 

HUfif + R^Sf^^-® 00 + L ±j 8 2 + R^] } , (38) 

where 8q CD and §^ USY -Q CD are the jjy divergent parts of the QCD and SUSY-QCD vertex 
corrections, respectively. These are given by 

S QCD = S SUSY-QCD = OC^Cf^ (gg) 

Air 

and 81,81, Sj 1 ^ and Sj 2 ^ are defined to be 

6, = (A Q R CD + Af S¥ ' QCD h mte , (40) 

8 2 = (Al CD + A S L US¥ - QCD )^ e , (41) 

= -R 2j - R 2j cos 29 - R 1A sin 29, (42) 

and 

Sf = L 2j - L 2j cos26 - L ljS m26. (43) 

We have checked analytically that all the ultraviolet divergences indeed cancel in the virtual 
corrections to the decay width, but the infrared divergent terms presist. 



4. Real corrections 

As is well known[13], to cancel the infrared divergences in the virtual corrections one 
needs to include real gluon emission, namely, t — > tix®g, as shown in Figs.l(e,f). As above, 
we will regulate the infrared divergences associated with the soft and collinear real gluon 
emission by the same finite small gluon mass A. In the calculation of the corrections due to 
real gluon emission to the partial width, it was necessary to perform the integration over the 



7 



three-body phase space. After tedious but straightfoward calculations we obtained 

STreai = "f/ 2 ^ {(l^f + \R^\ 2 )[I + % - 2(m* - (m|o - m|) 2 )/ i 
+2(m| - m 2 - m?o)(/ + A + m 2 / o + m |-^n)] 

+8m t m^oi2e(L* J -i2i iJ -)[(TO|o - m 2 - ra|)J i - - m 2 / o - io - A] j ,(44) 

Here we adopt the notation of Ref. [14] where the definition of the functions ij, Iij(m t , , tti^o) 
can be found. We also have checked numerically that the infrared divergences in <5r rea / and 
ST virt do indeed cancel. 



5. Numerical results and discussions 

In the following we give the numerical results for t — > tiXii where Xi is the lightest 
neutralino. In our numerical calculation we fixed M = 200 GeV, \x = —100 GeV and 
we used the relation M' = [7] to fix M' . For the parameters in stop sector we 

assumed Mf = M^ L and took the combination A t + fi cot f3 to be one parameter. Note 
that (A t + \i cot P) — corresponds to the case of no mixing in the stop mass matrix, 
Eq.(ll). There are then three free parameters in the stop sector and we chose m^, tan/3 
, and (A t + /xcot/3) as the three independent parameters. Other input parameters are 
m z = 91.188GeV, a em = 1/128.8, and G F = 1.166372 x 10' 5 (GeV)- 2 . The W mass was 
determined from [15] 

2 Q _ ™W\ _ ™ 1 U5) 

mw[1 Mf ~~ \/2Gf 1 — Ar ' (45) 
where, for a heavy top, Ar is given by [16] 

Figure 3 shows the relative correction to the decay rate ST/To, To being the tree-level 
rate, as a function of the lighter stop mass assuming trig = 500GeV and tan/3 = 11. The 
solid curve in Fig. 3 corresponds to A t + /i cot (5 — 0, the no-mixing case, while the dotted 
curve corresponds to A t + /i cot (3 = lOOGeV, a mixing case. Note that in Fig. 3 the lightest 
neutralino mass is m^o = 68 GeV. It is clear that the correction in the mixing case is 



larger than in the no-mixing case and can reach -20% for rrif = 100 GeV. Figure 4 shows 
the dependence of the relative correction to the decay width on the value of gluino mass 
for m~ tx = 50 GeV. Other parameter values are the same as in Fig.3. For the solid curve 
m~ tl = 50 GeV and m 4 ~ 2 = 64 GeV and there are two peaks at m g = 112 GeV and m g = 126 
GeV due to the fact that m t = 176 GeV and the threshod for open top decay into gluino 
and stop is crossed in these regions. For the dashed curve = 50 GeV and mf 2 = 194 
GeV and there is only one peak at m g = 126 GeV. When the gluino mass is heavier than 
200 GeV the correction in the mixing case is larger than in the no-mixing case and both 
corrections increase with gluino mass. Decoupling effects do not occur here, in contrast to 
the virtual SUSY corections to the decay and production processes in the SM. In Figure 5 
we present the dependence of the relative correction to the decay width on the value of tan (3 
assuming m g = 500GeV, = 50 GeV and A t + fi cot (3 = 100 GeV. Only in the region 
where tan f3 < 2 is the correction to the decay width very sensitive to the value of tan (3. 

In conclusion, we have shown that the one-loop QCD and SUSY-QCD corrections to 
t — > t\x Q j can exceed -10% of the tree level partial width in both the no-mixing and the 
mixing case of stop masses, and these corrections are not sensitive to tan/5 for tan/5 > 2. 

Note added: While preparing this manuscript the preprint of A.Djouadi, W.Hollik and 
C.Junger ( [hep-ph / 9605340]) appeared where the QCD correction to the process t — > tiXj is 
also calculated. But Eq.(14) of their original paper were not correct. Very recently, in their 
revised version this mistake has been corrected by them. We thank W.Majerotto for useful 
communication. 

This work was supported in part by the U.S. Department of Energy, Division of High 
Energy Physics, under Grant No. DE-FG02-91-ER4086. 
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Figure Captions 



Fig.l Feynman diagrams for the tree-level process t — > tiX® and the QCD corrections. 
Fig. 2 Feynman diagrams for the SUSY-QCD corrections. 

Fig. 3 The relative correction ST/Tq to the decay rate as a function of the lighter stop 
mass assuming m g = 500GeV and tan/? =11. The solid and dotted curves correspond to 
A t + ii cot f3 = (no mixing) and A t + //cot/3 = lOOGeV (mixing), respectively. 

Fig. 4 The relative correction <5r/r to the decay rate as a function of the gluino mass 
assuming = 50GeV and tan (3 = 11. The solid and dotted curves correspond to A t + 
yucot (3 = (no mixing) and A t + //cot (3 = lOOGeV (mixing), respectively. 

Fig. 5 The relative correction ST/Tq to the decay rate as a function of tan/3 assuming 
m- g = 500GeV, m~ tl = 50 GeV and A t + p cot (3 = lOOGeV. 
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